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Abstract Within the context of hierarchical scenarios of gravitational structure formation 
we describe how an evolving hierarchy of voids evolves on the basis of two 
processes, the void-in-void process and the void-in-cloud process. The related 
analytical formulation in terms of a two-barrier excursion problem leads to a 
self-similarly evolving peaked void size distribution. 

1. Introduction: Excursions 

Hierarchical scenarios of structure formation have been very succesfull in un- 
derstanding the formation histories of gravitationally bound virialized haloes. 
Particularly compelling has been the formulation of a formalism in which the 
collapse and virialization of overdense dark matter halos within the context of 
hierarchical clustering can be treated on a fully analytical basis. This approach 
was originally proposed by Press & Schechter (1974), which found a particu- 
larly useful and versatile formulation and modification in the the excursion set 
formalism (Bond et al. 1991). 

It is based on the assumption that for a structure to reach a particular nonlin- 
ear evolutionary stage, such as complete gravitational collapse, the sole condi- 
tion is that its linearly extrapolated primordial density should attain a certain 
value. The canonic example is that of a spherical tophat overdensity collapsing 
once it reaches the collapse barrier 5 C « 1.69. The successive contributions 
to the local density by perturbations on a (mass) resolution scale S m may be 
represented in terms of a density perturbation random walk, the cumulative of 
all density fluctuations at a resolution scale smaller than S m . By identifying 
the largest scale at which the density passes through the barrier 5 C it is possible 
(1) to infer at any cosmic epoch the mass spectrum of collapsed halos and (2) 
to reconstruct the merging history of each halo (see Fig. top right). 

In this study we demonstrate that also the formation and evolution of foam- 
like patterns as a result of the gravitational growth of primordial density pertur- 
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Figure 1. Illustration of the two essential "void hierarchy modes": (top) the void-in-void 
process (top), with a void growing through the merging of two or more subvoids; (bottom) 
the void-in-cloud process: a void demolished through the gravitational collapse of embedding 
region. 

bations is liable to a succesfull description by the excursion set analysis. This 
is accomplished by resorting to a complementary view of clustering evolution 
in which we focus on the evolution of the voids in the Megaparsec galaxy and 
matter distribution, spatially the dominant component. An extensive descrip- 
tion of this work can be found in Sheth & Van de Weygaert (2003, MNRAS, 
submitted). 



2. Void Evolution: the two processes 

Primordial underdensities are the progenitors of voids. Because underdensi- 
ties are regions of suppressed gravitational attraction, representing an effective 
repulsive gravity, matter flows out of their interior and moves outward to the 
edges of these expanding voids. Voids expand, become increasingly empty and 
develop an increasingly spherical shape (Icke 1984). Matter from the void's 
interior piles up near the edge: usually a ridge forms around the void's rim and 
at a characteristic moment the void's interior shells take over the outer ones. 
At this shell-crossing epoch the void reaches maturity and becomes a nonlinear 
object expanding self-similarly, the implication being that the majority of ob- 
served voids is at or near this stage (Blumenthal et al. 1992). As voids develop 
from underdensities in the primordial cosmic density field, the interaction with 
internal substructure and external surrounding structures translates into a con- 
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tinuing process of hierarchical void evolution (Dubinski et al. 1993, Van de 
Weygaert & Van Kampen 1993). 

The evolution of voids resembles that of dark halos in that large voids form 
from mergers of smaller voids that have matured at an earlier cosmic time 
(Fig-E top row). However, in contrast to dark halos, the fate of voids is ruled 
by two processes. Crucial is the realization that the evolving void hierarchy 
does not only involve the void-in-void process but also an additional aspect, 
the void-in-cloud process. Small voids may not only merge into larger en- 
compassing underdensities, they may also disappear through collapse when 
embedded within a larger scale overdensity (Fig. ^ bottom row). In terms of 
the excursion set approach, it means that the one-barrier problem for the halo 
population has to be extended to a more complex two-barrier problem. Voids 
not only should ascertain themselves of having decreased their density below 
the void barrier S v of the shell-crossing transition. For their survival and sheer 
existence it is crucial that they take into account whether they are not situated 
within a collapsing overdensity on a larger scale which crossed the collapse 
barrier S c . They should follow a random walk path like type "3" in Fig.|2](bot- 
tom right), rather than the void-in-cloud path "4". The repercussions of this 
are far-reaching and it leads to a major modification of the void properties and 
distribution. 
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Figure 2. Left: Void Size Probability Function for various values of the void-and-cloud 
parameter S C /8 V : 1/2 (short dashed), 1 (dashed), 2 (dot-dashed), and oo (solid). For reference, 
the dotted curves are predictions by the "adaptive peaks model" of Appel & Jones (1990). 
Right: Excursion Set Random Walks, local density <5(x) as function of mass resolution S m - 
Clearly, the fluctuations are larger as resolution S m increases. (Top) Halo formation. Horizontal 
dashed bar is the collapse barrier S c . (Bottom) Fate of void V2, illustrated by 2 possible random 
walks: curve "4" is a void-in-cloud process, curve "3" a void-in-void process. 
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3. Void Distribution: Universal, Peaked and Self-Similar 

Analytically, the resulting expression follows by evaluating the fraction of 
walks which first cross S v at S, and which do not cross 5 C until after they have 
crossed 5 V (Sheth & Van de Weygaert 2003). An insightful approximation of 
this distribution, in terms of the "self-similar" void density u = 51/ S is 

in which the D = (\5 V \/(5 C — 5 V ) parameterizes the the relative impact of void 
and halo evolution on the hierarchically evolving population of voids. 

The resulting distributions, for various values of D, are shown in Figure |2] 
(left). The void size distribution is clearly peaked about a characteristic size: 
the void-in-cloud mechanism is responsible for the demise of a sizeable popu- 
lation of small voids. The halo mass distribution diverges towards small scale 
masses, so that in terms of numbers the halo population is dominated by small 
mass objects. The void population, on the other hand, is "void" of small voids 
and has a sharp small-scale cut-off. 

Four additional major observations readily follow from this analysis: (•) 
At any one cosmic epoch we may identify a characteristic void size which 
also explains why in the present-day foamlike spatial galaxy distribution voids 
of ~ 20 — 30/i _1 Mpc are the predominant feature; (•) The void distribution 
evolves self-similarly and the characteristic void size increases with time: the 
larger voids present at late times formed from mergers of smaller voids which 
constituted the dominating features at earlier epochs (Fig.^ top panels); (•) 
Volume integration shows that at any given time the population of voids ap- 
proximately fill space, apparently squeezing the migrating high-density matter 
in between them; (•) As the size of the major share of voids will be in the 
order of that of the characteristic void size this observation implies that the 
cosmic matter distribution resembles a foamlike packing of spherical voids of 
approximately similar size and excess expansion rate. 

In all, a slight extension and elaboration on the original extension formulism 
enables the framing of an analytical theory explaining how the characteris- 
tic observed weblike Megaparsec scale galaxy distribution, and the equiva- 
lent frothy spatial matter distribution seen to form in computer simulations of 
cosmic structure formation, are natural products of a hierarchical process of 
gravitational clustering. A continuously evolving hierarchy of voids produces 
a dynamical foamlike pattern whose characteristic dimension grows continu- 
ously along with the evolution of cosmic structure, a Universe which at any 
one cosmic epoch is filled with bubbles whose size corresponds to the scale 
just reaching maturity. 



